^^^wmm^m 


o    'It 


'.^S;> 


'C^ 


V  ■ 


■*:  -' 


hi"  y 


^  ■'!'J^*' 


y^w 


»". 


LIBRARY  OF  THE 

UNIVERSITY  OF  ILLINOIS 

AT  URBANA-CHAMPAIGN 

510.84 

no. 100-110 
cop. 3 


M 


ho. '05 
Cop .  3 


UNIVERSITY  OF  ILLINOIS 
GRADUATE  COLLEGE 
DIGITAL  COMPUTER  LABORATORY 


REPORT  NO,  105 
A  CRITICAL  STUDY  OF  AlWL 
by 
Herman  Bottenbruch 


December  8,  I96O 


This  work  was  supported  in  part  by  the 
Atomic  Energy  Commission  and  the  Office  of  Naval  Research 
under  Contract  AT(ll-l)-4l5 


^^e-^  FOREWORD 


In  February ;,  I96O,  Dr,  Herman  Bottenbruch  from  Damstadt,  Germany* 
gave  a  seminar  on  ALGOL  at  the  Digital  Computer  Laboratory.  This  seminar  was 
followed  by  an  active  discussion,  in  which  Dr.  Bottenbruch  and  various  members 
of  the  Laboratory  participated.  Some  of  the  points  raised  in  this  discussion 
made  it  obvious  that  it  would  be  worthwhile  to  examine  ALGOL  in  more  detail, 
with  special  attention  to  properly  evaluating  its  merits  and  its  deficiencies. 
Such  an  effort  was  warranted  to  be  especially  valuable  due  to  the  fact  that 
most  of  the  published  literature  on  ALGOL  does  not  examine  these  questions 
adequately  or  else  treats  them  in  an  unscientific  or  emotional  fashion. 

Dr.  Bottenbruch  kindly  spent  the  first  two  weeks  of  July  in  the 
Digital  Computer  Laboratory,   During  this  period,  he  presented  more  detailed 
descriptions  of  ALGOL  and  entered  into  discussions  concerning  the  general  ques- 
tions of  ALGOL  programming  procedures  and  the  machine  independence  of  ALGOL. 
He  prepared  some  ALGOL  programs  for  some  of  the  processes  suggested  by  members 
of  the  Laboratory  as  examples  illustrating  the  aforementioned  questions. 
Some  of  these  are  contained  in  this  report. 

This  report  is  a  distillation  of  those  discussions.  The  Laboratory 
assumes  no  responsibility  for  the  opinions  and  ideas  expressed  there  and 
recorded  here.  These  have  been  influenced  by  the  criticisms  and  suggestions 
which  were  made- -hopefully  in  the  right  direction. 


*This  year  visiting  with  the  Mathematics  Panel  of  ORNL,  Oak  Ridge,  Tennessee, 
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I.   GEEERAL  REMARKS  ON  PROGRAMMING  DIGITAL  COMPUTERS 


1,1  Flow  Diagram^  Machine  Program  and  Algorithm 

In  the  very  first  paper  on  programming  a  general  purpose  digital  computer 
[l]  by  Goldstine  and  von  Neumann,  the  importance  of  making  a  "plan"  of  the  computation 
before  writing  the  detailed  ma,chine  program  in  ma,chine  code  was  stressed.  In  their 
philosophy,,  the  construction  of  the  program  consisted  of  the  following  steps: 

Step  1.  General  layout  of  the  process  to  be. carried  out,  written  in 
ordinary  English  with  the  use  of  mathematical  formulae. 
This  general  plan  also  made  reference  to  storage  locations 
in  symbolic  form. 

Step  2.  Rigid  statement  of  the  process  to  be  carried  out.   In  [l]  the 

"flow-diagram"  is  defined  to  be  used  as  the  language  for  stating 
this  second  phase  of  programming. 

Step  3»  Translation  of  the  flow-diagram  into  machine  code.  This  third 
step  was  accomplished  in  two  sub-steps: 

Step  3a«  Writing  machine  code  in  symbolic  addresses. 
Step  3t).  Writing  machine  code  in  absolute  addresses. 

Both  the  notation  used  in  the  flow-diagram  of  [l]  as  well  as  the  way  in  which  the 

flow-diagrams  for  particular  problems  were  constructed,  were  essentially  machine 

1  2 
independent.  '        The  authors  of  [l]  expressed  in  several  places  their  opinion 

that  the  essential  steps  of  programming  planning  were  completed  after  Step  2,  and 

that  the  machine  coding  was  a  "straight  forward"  and  "simple"  task.  They  did, 

however,  not  bother  about  a  computing  machine  to  do  this  simple  job. 

The  second  major  contribution  to  the  art  of  programming  [2]  expressed  - 
not  explicitly,  but  in  practice  -  a  rather  different  point  of  view.  The  authors 
held,  that  Step  2,  the  rigid  statement  of  the  procedure,  should  as  far  as  possible 
take  advantage  of  the  order  code  of  the  particular  machine  to  execute  the  procedure. 
So  for  them  the  formulation  of  the  program  and  translation  into  machine  code  collapsed 


1.  An  exception  to  this  is  the  fact,  that  the  "alternative -boxes"  made  the  choice 
of  th.e  continuation  of  path  dependent  on  the  sign  of  an  algebraic  expression 
rather  than  the  value  of  a  boolean  expression, 

2.  See  however,  the  remarks  on  "machine  independence"  in  1.3. 
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into  one  step;  Formulation  was  done  in  machine  code.  Tbey   devised,  however, 
several  methods  to  make  this  "formulation"  easier  for  h-i;aT..an  beings,  such  as  relative 
addi'esses,  preset  parameters,  automatic  memory  assignment  to  sub-routines  etc.  They 
were  indeed  the  first  to  stress  the  poin^  that  the  program  fed  into  the  machine  and 
the  program  to  be  executed  by  the  machine  need  not  be  exactly  alike  -  the  machine 
is  able  to  perform  certain  elementary  transformations  on  the  external  program 
fed  into  it. 

The  third  approach  to  programming  was  given  by  Rutisha,user  [3],  following 
early  ideas  of  Zuse  [h],   which  in  certain  respects  seems  to  be  a  revival  of  the  von 
Ne^jmann  approach  and  in  certain  o^her  respects  seems  a  natural  extension  of  the  ideas 
of  Wilkes,  et.al.  First  it  stresses  the  point  tnat  before  the  actual  coding  begins  a 
rigid  plan  of  the  computation  (Rechenplaa)  has  to  be  constructed.  But  the  language  in 
which  to  foi'^ulate  this  plan  is  chosen  such  that  a  general  purpose  computer  can  be 
programmed  to  translate  that  plaxi  into  its  own,  order  code.  The  obvious  advantage 
which  cne   gets  in  having  a  language  which  can  be  translated  by  a  machine  is  in  part 
offset  by  the  fact  that  some  important  characteristics  of  the  von  Netunarin  flow- 
diagram  are  lost;  e.g.  the  "bound  variables",  ar.d  the  "storage  tables"  with  their 
separation  of  mathematical  quantities  and  storage  locations.  On   the  other  side,  the 
obvious  advantage  which  one  gets  in  having  a  programming  language  which  allows  one  to 
write  programs  which  look  more  like  computa,tion  plans  thaxi  machine  programs  is  in 
part  offset  by  the  fact  the  progrsaiimer  no  longer  has  full  control  over  the 
facilities  of  a  particular  machine  -  and  It  is  of  co;^^Be  understood  that  no  trans - 
la,tlon  program  will  do  as  good  a  Job  of  exploiting  the  possibilities  of  a  computer 
as  an  experienced  progra;mner.  The  ideas  of  Laning  and  Zierler  are  along  the  same 
line,  and  with  FORTRAN  [6]  this  method  became  a  m.ost  powerful  tool  in  programming. 
ALGOL  [7],  [8]  is  the  attempt  of  many  of  the  groups  interested  in  this  approach  to 
programming  to  establish  a  laniversal  language  suitable  for  translation  into  the 
order  code  of  any  general  purpose  computer. 

1.2  Is  ALGOL  Ma,chine  Independent  and  Universal? 

ALGOL  is  designed  to  describe  computation  processes  on  present  day 
computers.  ALGOL  programs  should  be  "understood"  by  human  beings  as  well  as  by 
machines.  In  this  paragraph  we  describe  the  potentialities  of  ALGOL  of  describing 
computation  processes  more  or  less  m  abstracto  -  without  particularly  referring 
to  men  or  machines  to  understand  these  descrlp-rions.  In  1.3  we  discuss  how  good 
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ALGOL  is  for  understanding  by  a  hiunan  teing^  and  in  l.k   we  discuss  how  good  ALGOL 
is  for  decipliering  'by  a  machine.  Although  all  three  questions  -  adequateness  for 
expression  ,    understanding  by  humans  and  deciphering  by  machines  are  closely 
related;,  a  separate  discussion  is  at  least  easier  for  the  writer  of  this  report. 

Since  a  general  purpose  computer  can  simulate  any  other  general  purpose 
computer;  the  order  code  of  any  general  purpose  computer  is  in  a  certain  sense 
"machine  independent",  ALGOL  is  called  machine  independent  in  a  different  sense. 
To  explain  this  more  fully  let  us  distinguish  between  "what  is  done"  on  computing 
machines  and  "how  this  is  expressed". 

Now  the  "what  is  done"  on  different  computers  is  never  the  same  if  it  is 
described  in  tenns  of  the  order  code  for  those  computers.  If  however  it  is  described 
only  in  more  general  terms  like  "a  loop  is  formed"  or  "an  address  is  computed",  then 
in  these  general  terms  the  "what  is  done"  may  well  be  the  same,  and  a  universal 
language  on  some  level  is  siirely  possible.  The  fort 'Xiate  fact  about  ALGOL  is  that  it 
is  at  the  same  time  general  enough  that  using  its  terms  the  description  of  "what  is 
done"  is  identical  for  different  computers,  but  on  the  other  hand  is  close  enough 
to  "machine  thinking"  that  translation  into  the  order  code  of  computers  is  not  too 
difficult  and  yet  gives  efficient  codes.  So  If  the  statement  "ALGOL  is  machine 
independexit"  means  anything  it  means  this;   In  terms  of  ALGOL  the  descriptions  of 
a  particular  computing  process  for  different  machines  are  identical,  even  comparing 
the  most  efficient  formulations  of  the  procedure  for  each  computer. 

Now  this  statement  is  not  a  mathematical  theorem,  and  its  validity  lies 
somewhere  between  true  and  false,  depending  upon  the  class  of  machine  and  computation 
procedures  under  consideration. 

Let  us  consider  several  examples,  in  which  ALGOL  proves  to  be  more  or  less 
machine  independent . 


Example  Is 


boolean  array  a,  b  [itkO]^ 
for  ii      =  1  step  1  until  ^0  do 
a[i3s=  a[i]  -  b[i],° 


3,  This  is  different  from  saying  that  the  order  code  for  a  particular  universal 
computer  is  machine  independent. 
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This  may  be  a  perfect  description  of  a  computing  process  for  many  machines, 
especially  decimal  machines.  Many  machines,  however,  can  store  the  values  of  the 
arrays  a  and  b  in  one  (binary)  computer  word  and  in  one  single  instruction  can  form 
i-he  logical  sum  of  two  such  words  as  required  in  the  program.  For  those  machines, 
apparently,  the  ALGOL  program  is  not  good. 

Example  2: 

10 

E   a,  b 

i-1 

cos  a  =   ,     Hyr 


1   1 


yczap  •  (Eb^) 


ALGOL  Program  Is 

SI;  =  S2;  =  S3i   =   0; 

for  i ;  =  1  step  1  until  10  do  begin 

SI;  =  SI  +  a[i]  x  b[l]5 

S2;  =  S2  +  a[i]  x  a[i]; 

S3;  =  S3  +  b[l]  X  b[i]  end; 
cos;  =  Sl/sqrt  (S2  x  S3)j 

ALGOL  Program  2: 

SI;  =  S2;  =  S3;  =  0; 

for  i;  =  1  step  1  until  10  do  SI;  =  SI  ■♦-  a[i]  x  b[ij; 
for  i;  =  1  step  1  until  10  do  S2:  =  S2  +  a[i]  x  a[i]j 
for  i;  =  1  step  1  until  10  do  S3;  =  S3  +  b[i]  x  b[i]j 
cos:  =  Sl/sqrt  (S2  x  S3)i 

Programs  1  and  2  are  functionally  equivalent)  they  only  operate  different,, 
some  machines  formulation  1  might  be  preferable,  on  other  machines,  which  have  a 
sma.ll  high-speed  memory  that  can  accomodate  a  small  loop,  formulation  2  might  be 
preferable.  This  seems  to  be  an  example  that  shows;  Even  if  the  "what  is  done" 
is  not  the  same  for  different  machines,  it  still  is  of  great  advantage  to  have  only 
one  language  describing  the  two  different  possibilities  -  if  that  language  is 
universal . 


The  machine  independence  claim  could  be  maintained  by  saying  that  a  translation 
process  is  possible  which  translates  that  ALGOL  program  into  the  described 
machine  program.  This  woijild  make  the  questions  under  discussion  depend  on  the 
existence  or  possible  existence  of  translation  processes.  The  role  of  trans- 
lators is  more  fully  discussed  in  l^k-. 
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There  might  be  ma.chlnes  with  three  different  computing  units  forming  the 
three  sims  in  parallel.  For  a,  description  of  this  procedure  ALGOL  is  not  clearly 
good  heca.use  there  is  no  means  in  it  to  describe  parallelism. 

So  from  the  continuation  of  possible  cases  we  can  point  to  the  following 
a.s  being  essentially  different: 

1)  An  ALGOL  program  may  be  a,  good  description  for  computing 
procedures  on  machines  A  and  B. 

2)  An  ALGOL  program  which  is  good  for  machine  A  may  be  poor 
for  machine  B. 

3)  A  procedure  for  ma,chine  A  cannot  be  described  in  ALGOL. 

The  question  of  whether  ALGOL  is  "universal"  or  not  is  of  a  different 
kind  and  cannot  be  answered  by  logical  reasoning  but  only  by  observation  of 
future  events .   Certainly  even  if  it  were  completely  ma.chine  independent  but 
were  not  universally  a.ccepted  it  would  not  be  universal^  and  if  a.  language  of 
a  particularly  cleverly  constructed  computer  would  be  universally  a.ccepted  it 
would  be  universal  -  even  though  "machine  independent".   So  logically  the  two 
properties  do  not  seem  to  be  connected.  Actually^  they  are^,  however,  because 
only  a  language  which  is  rea.sonably  machine  independent  will  have  a,  chance  to 
be  universally  a.ccepted  for  the  description  of  computing  processes. 

1.3  ALGOL,  Common  La.ngua,ge  and  Mathematical  Nota,tion 

The  fundamental  difference  between  ALGOL  and  ma.thema.tica.l  nota.tion  is 

that  ALGOL  is  dynamic;  (all  the  statements  of  ALGOL  are  of  the  kind:   Do  something) 

a.nd  mathema.tics  is  sta.tic  (assertions  a.re  ma.de).   A  ma.thema.tical  sta.tement  is  e.g. 

2x2  =  4,  which  certainly  does  not  imply  tha,t  something  must  be  computed.  The 

sta.tement!   if  X.  .,  =77  (zr-   +  X.)  for  all  i.  then  X.  converges  towards  ,/a  does 
1+1   2  ^X.    1  '       1       ^  -^ 

not  imply  that  something  has  to  be  computed,  but  states  that  a  certain  rela.tion 
holds  (convergence  of  a.  sequence)  if  certa.in  other  relations  hold  (between  the 
elements  of  the  sequence).   Common  language  has  to  be  used  in  addition  to  ma.the- 
ma.tical notation  to  sta,te  computation  processes,   e.g.:   "Set  X.  =  a,  and  compute 

1    a  "*" 
X.  -,  from  X.  for  several  i  a.ccording  to  the  rule  X.  ,  =  7:  x  (zr-   +  X.)»   Continue 
1+1       1  ^ 1+1   2   ^Xj_    1' 

until  two  successive  values  X.  ,  and  X.  are  sufficiently  close  together.   Use 
only  one  storage  location  for  the  different  X.." 


5.  A  limited  amount  of  parallelism  is  in  arithmetic  expressions.   For  the  expression 
(a[i]  +  b[i])  X  (c[j]  +  d[j]),  there  is  nothing  tha.t  states  which  of  the  four 
a.ddresses  involved  has  to  be  computed  first,  so  they  ma.y  be  computed  in  pa.ra.llel. 
Also,  the  two  sums  which  form  the  factors  ma.y  be  computed  in  pa.rallel. 
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The  corresponding  ALGOL  statement  would  read: 

X:   =  aj 

L:   d:   =  0.5  x  (|  +  X); 

if  abs  (d^l*  -)  >  2  X  (_   then  go  to  L; 

comments  at  this  point  ahs  {sqrt   (a)  -  X)  <  x.} 

Note  that  in  the  ALGOL  program  no  reference  is  made  to  the  "mathematical"  quantities 

1^  2'      S"*  *  *  *  * 

Apparently  the  ALGOL  program  is  more  precise,  but  still  there  is  -  for 
people  used  to  reading  ALGOL  -  no  price  paid  for  this  precision  as  far  as  under- 
standing the  program  is  concerned. 

We  will  discuss  the  question;  Does  the  precision  of  ALGOL  programs  and 
the  restriction  Imposed  on  the  use  of  terms  affect  the  understanding  of  computation 
processes  described  in  ALGOL? 

Let  us  consider  the  following  ALGOL  procedure  for  numerical  integration 
by  Simpson's  rule  with  automatically  adjustable  step  length  (compare  ["j])- 

procedure  Simpson  (a,b,f,  Qj;  value  a,b,0  real  a,b,£; 

begin  n;  =  Ij  h:  =  (b-a)/3| 

I:  =  (b-a)/6  X  (f(a)  +  k   x  f(^)  +  f(b))| 
St   =   f(a)  +  2  X  f((a+b)/2)  +  f(b); 
Lt      n;  =  2  X  nj  h:  =  h/2| 

T:  =  0,"  for  i;  =  1  step  2  until  2  x  n  do  T:  =  T-i-f(a+l  x  h); 
NI;  =  (S+1^  X  T)  X  h/3i 

if  abs  (NI  -  I)  >  £  then  begin  I:  =  NI;  S:  =  S  +  2  x  T;  go  to  L  end; 
Simpson:  =  NI  end. 

It  is  not  altogether  obvious,  even  for  an  experienced  ALGOL  programmer, 
to  figure  out  how  that  program  achieves  its  aims  and  whether  it  does  achieve  those. 
The  following  verbal  description  seems  preferable  as  far  as  understanding  is 
concerned 2 

/^^ 

/     f(X)  dX 

—^  a 

is  formed  using  Simpson's  rule  by  dividing  the  Interval  from  a  to  b  into  1,2,^,8, .. .,2' 

sxeps.  Let  the  n   approximation  be  called  I  |  then  the  process  is  stopped  if 

-  I  -,  I  <  f  .   In  each  step  the  quantities  S  and  T  are  formed,  where: 
n    n-1'    ^  jr-     -u  n     n  ' 
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J. 


n 


i-1,3. 


f(s.  +  -- 
2 


i) 


ilL       "X 


:a)  +  f(l5)  +  2  X    Z 


fCa  + 


"b-a 


2 


n 


i) 


T     is  fonrmed  by  actually  eral-oating  the  sibils  vliereas  S     is  computed  as; 
n  n 


n 


5     .   +  2  :x  T   , 
ii-l  n 


Fr>om  S     and  T     we  get  I     "by; 
n  n        *^         n    "^ 


I     =.  (S     +  2  x  T   )    •  ^^ 


"b— a 
The  q'aantities  —  are  of  course  computed  recuTBiTelVo     TKiere  is  only  one  storage 

location  for  all  the  S  3  aao'b'.h.ejT  one  for  all  t'lie  T  , 


another  for  all  the  I 


By  this  procedur'e,  we  evaloate  the  fanietioa:.  f  oaly  once  for  every  point 
used  in  the  f.lnaJL  app.ToxiLiEationj,  wheireas  we  wo'aJd  ewBl-a8.te  the  fiaaction  more  often 
in  a  more  straight  forwstifd  proced-ore^ 

The  constractioa  of  an  ALGOL  prog'jr^js.  acc!^rd3:«ig  to  the  given  description 
is  simple  -    maybe   even  simpler  tibiati  uaderetaadijfflg  the  ALSOL  prograia  for  the  same 

procedu^'e  given  above.     There  a.re  other  :lnstaaces^  where  the  verb&l  description 
of  a  computation  process  1b  of  almost  uo  laelp  ia  Vrit:lng  the  corresponding  ALGOL 
program o     In  these  cases j,   coanEOja  lau^aage  is  a.ot  adeqmte  for  eacpressing  the 
details  of  the  computation  wheresjs  ALGOL  still  is.     Eramples  are  the  matrix- 

sortiag  prog.r?uii  and  the  multiplication  of  two  sparse  matrices  in  Chapter  2o 


Here,  ALGOL  -  or 


Mug,   eejiiiTalent  »  is  the  leas^Qiage  to  describe  computations, 


AlsOj,  in  these  cases,  of  nourse^,  a  verbal  statement  of  the  program  is  of 
help  in  understanding  iChe  ALGOL  program.  The  art  of  ili^jBLinating  compu- 
tation processes^,  whether  in  ALGOL  or  in  mackine  code,  is^  however^  in  a 
very  poor  condition, 
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1.4  ALGOL  arid  Translators 

The   question  of  whether  ALGOL  is  adequate  to  describe  efficient 
co3ipui'.at.ion  processes  for  computers  depends  on  the  translation  process.  If 
translation.  i.&   doue  l.y  irran,  aolsody  can  discuss  the  question  because  nobody  knows 
how  somebody  is  going  to  translate  a  particular  program.  So  we  restrict  dis- 
cissi.on  to  translation  processes  which  can  be  carried  out  by  machines.  However, 
since  aay  translation  process  can  be  programmed  for  a  computer  we  restrict  dis- 
CuSBlon  even  further  to  "reasonable  translators",  that  is  those  which  can  be 
constraetad  without  undue  effort,  which  work  fast,  or  -  in  short  -  which  exist 
or  will  come  fortia  In.  the  near  future.  This  is  still  very  vague.  It  is  however, 
prooat-ly  tr^^e  tMt  tlie  border  line  for  a  reasonable  translator  lies  in  the  tran- 
sition from  a  "static"  understanding  to  the  "dynamic"  understanding  of  a  program. 
So  e.g.  a  "reasonable"  trenslator  might  be  able  to  find  out  that  one  of  the  factors 
inf clTsd  in  a  multiplication  is  a  constant  and  a  power  of  two,  and  it  might  code 
tMs.  ffi-iltipli cation  as  a  shift  ratiier  than  an  actual  multiplication,  if  shifting  is 
fast  and  miiltiplisatioa  is  ?iOt,  A  reasonable  translator  might  not  be  able  to  find 
ail  insportant  fact  In  the  fol.iDwi;ng  ALoOL  programs 

a[l]s  -  1| 
for  is  =^  2  step  1  .ant  11  50  do 
a[i];  =^-  i  +  a[i-l]> 

The  impcroaat  fa  at  is,  that  the  result  a[:i3  is  used  as  one  of  the  operands  in  the 
ii.ext  pas&s  thro'ugSii  the  loop.  So  an  intelligent  taranslator  (in  fact,  any  himian  being) 
woTild  lea^B  the  result  in  -nhe  s,c.rj!LjilB,to7j   t.o  haye  it  available  next  time.  A 

reasoiaable  treaislator  woiild  not  fee  able  to  do  this. 

Of  course,  there  is  the  possibility  of  expressing  the  dynamic  fact  "a[i] 

is  'ased  neiK"*'.  time"  by  a  tiaich  of  static  i^-zl.es  aboiit  sequences  of  symbols  which 
represent  p:i?og:rsflis  for  which  tJjis  dynamic  fact  is  tr^ie.  It  is  however,  not  simple 
to  "t'kJak  out"  the  rTil.es,  and  the  coding  of  these  rules  for  computers  is  even 

harder  o 

Another  eifSUEple  of  "dynamic  understanding"  necessary  for  optimisation 
is  fixe  eaeaaLple  in  1,2 »  Here  it  must  be  fo^aid  out  that  ail  the  kO  executions  of 
the  loop  can  I'.e  .^s.^rrted  c  :,t  ijo.  parallel  ->  a  fact  which  is  generally  not  true  for 
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There  are  a  few  instances  where  dynamic  changes  in  an  ALGOL  program 
during  translation  are  almost  necessary.  These  are? 

1)  Elimination  of  common  expressions,  and 

2)  Recursive  computation  of  simple  expressions. 

To  l)t     In  a  statement 

X;  =  (a  X  fc  +i)  /  (a  X  b  +  g) 

it  caui  be  found  out  using  rather  simple  rules  that  the  values  of  a  and  b 
do  not  change  from  the  first  evaluation  of  the  product  to  the  second.  So 
a  translator  would  produce  a  machine  program  that  compxites  a  x  b  only  once 
and  stores  it  for  fiirther  reference.  If  in  a  sequence  of  statements,  none 
of  which  is  labeled^  the  product  a  x  b  occurs  in  several  places,  and  if  a 
and  b  do  not  occur  on  the  left  side  of  an  assignment  statement  in  the  given 
sequence,  then  the  product  may  be  precamputed.  If  a  product  a  x  b  occurs  in 
a  loop,  and  if  neither  a  nor  b  are  changed  in  the  loop,  then  a  x  b  may  be 
computed  outside  the  loop. 

There  are  more  complicated  cases,  in  which  an  expression  occurring  in 
several  places  may  be  precomputed,  which  are  not  covered  by  the  foregoing 
discussion. 

Consider  the  loop 

for  j;  =  1  step  1  until  10  d£ 
a[j];  =  10  X  j  +  ii-i 

The  computation  of  the  product  10  x  j  by  actual  multiplication  is 
rather  foolish.  The  following  program  is  betters 

a[l]:  =  l4.j 
for  j;  =  2  step  1  until  10  do 

a[j];  =  a[j-l]  +  lOj 

What  does  "better"  mean  exactly?  If  the  ALGOL  translator  for  a  machine 
cstn  find  out  that  the  q^mntities  a[j]  can  be  formed  recursively  using 
only  additions,  it  would  be  immaterial  which  of  the  above  is  written. 

To  find  out  that  several  quantities  can  be  computed  recursively 
clearly  involves  some  kind  of  dynamic  understanding,  and  a  reasonable 
translator  should  have  some  provisions  tx)  find  the  most  important  schemes 
for  recursive  calculation.  It  must  be  clear  that  the  incorporation  of 
optimisation  features  into  translators  may  be  of  importance  for  the 
judgment  of  an  ALGOL  program. 
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lo5  Philosophy  and  Politics 

It  is  felt  by  many  people  that  present  life  with  its  many  specializations 
is  complicated  enough.  While  further  specialization  and  division  of  labor  is  still 
efficient,  it  has  long  ceased  to  he  pleasant,  and  some  people  think  it  is  even 
worth  some  money  to  prevent  further  specialization  with  its  automatic  construction 
of  walls  between  closely  related  fields.  Some  even  hold  that  this  -  in  the  long 
run  -  will  pay  off  in  dollars  or  equivalent. 

The  advent  of  digital  computers  had  two  major  effects:  It  saved  money 
at  many  places  where  these  giants  do  essentially  the  same  job  as  men,  but  do  it 
cheaperi  and  it  pushes  the  border  line  of  human  activities  into  new  regions  thus 
creating  new  challanges  to  mankind.  But,  alas,  to  get  a  machine  to  do  what  one 
wants  it  to  do  is  a  formidable  job,  even  if  one  is  (or  believes  to  be)  perfectly 
sure  about  what  one  wants  it  to  do.  Furthermore,  the  situation  gets  worse  by  the 
fact  that  the  same  problem,  done  in  essentially  the  same  way  on  two  different 
machines  may  look  quite  different,  so  that  one  machine  can  hardly  take  advantage 
of  codes  written  for  another.  This  situtation  then  automatically  creates  this  new 
specialization,  the  programmer,  intending  to  bridge  the  ga,p  between  the  man  with  a 
problem  and  the  machine,  sometimes  proves  to  be  more  of  a  wall  between  these  two. 

A  suggestion  to  remedy  this  situtation  is  to  create  a  language,  which: 

1)  States  only  the  essentials  of  a  computing  procedure,  leaving 
aside  the  details  necessary  for  a  machine  code, 

2)  Includes  all  the  details  which  are  necessary  to  state  a 
coDQiutatioii  procedure  in  a  rigid  way,  thus  making  it 
possible  to  translate  it  mechanically  into  machine  program, 

3)  Because  l)  and  2)  may  become  universally  accepted  as  a 
vehicle  for  stating  computation  processes,  this  provides 
the  advantage  of j 

a.  taking  the  drudgery  out  of  programming, 

b.  making  programs  written  for  one  machine 
accessible  to  another  machine,  and 

c.  removing  walls  between  man  and  machine 
and  man  and  man. 
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One  partial  realization  of  this  suggestion  is  ALGO|jo  Now  it  must  be 
proved  that  the  idea  is  practical »  The  answer  rests  very  much  on  the  fact  of 
whether  it  really  can  state  the  "essentials"  of  a,  computation  procedure.   If 
some  little  detail  may  speed  up  com.putation  by  a  factor  of  5  or  10^,  that  "detail" 
surely  is  essential,, 

The  question  of  essentials  is  discussed  more  fully  in  the  other  parts 
of  this  report c  In  some  cases,  it  turns  out  that  AKSOL  cannot  adequately  describe 
some  details  which  on  some  machines  may  be  essential  as  far  as  computing  time  is 
concerned o  This  is  just  too  bad  -  whether  for  ALGOL  or  for  those  machines  is  up 
to  the  future. 

Sbme  price  for  leaving  out  the  formulation  of  a  computation  process 
must  be  paid  on  any  asachine.  It  is  not  designed  to  meet  specifically  the  needs 
of  ALGOL^  which  after  all  might  not  be  a  bad  idea  after  things  have  become  more 
settled.  If  it  were  only  for  the  philosophical  satisfaction  of  keeping  the  world 
from  getting  more  complicated  not  many  people  would  be  willing  to  pay  this  price. 
There  axe,  however,  some  savings » 

Pirst^  program  preparation  cost  money^  and  simplified  program 
preparation  together  with  program  exchange  will  reduce  these  costs.  Hopefully, 
too,  the  better  programs  (those  which  ajre  cheaper  in  execution)  will  have  a 
higher  transfer  rate,  thus  saving  money  and  effort  also  on  this   side. 
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II.  THE  SCOPE  OF  ALGOL 

This  report  assimes  that  the  reader  knows  the  ALGOL  language,  or  at 
least  is  able  to  understand  a  program  written  in  ALGOL.  Section  2,X   contains 
a  short  review  of  ALGiOL,  which  may  have  the  advantage  over  the  official  report 
of  being  more  readable  -  at  the  price  of  not  answering  all  questions.  In 
Sections  2.2  and  2.3  some  ALGOL  programs  are  given.  While  everybody  believes 
that  ALGOL  is  good  enough  for  processes  which  consist  mainly  of  straightforward 
arithmetic  calculation^,  the  adequateness  of  ALGOL  for  stating  procedures  which 
consist  mainly  in  computation  of  addresses  and  moving  numbers  around  is  some- 
times questioned o  Therefore^,  the  examples  in  2.2  and  2.3  have  been  chosen  from 
this  latter  class. 

2.1  Short  Review  of  ALGOL  60 

In  an  ALGOL  program  a  programmer  can  refer  to  a.  storage  location  in 
symbolic  form.  A  storage  location  may  be  referred  to  by  a  single  name  (identifier) 
for  only  this  storage  location,  or  it  may  be  the  element  of  an  array.   In  the 
la,tter  case,  the  array  has  one  single  name,  and  the  individual  elements  of  the 
array  are  distinguished  by  subscripting.  The  size  of  the  array  must  be  indicated 
in  the  ALGOL  program  by  an  array  declaration  (See  below)  which  e.g.  may  have  the 
form:   array  a.[l:10],  b[l;4,  3::5]>  with  the  meaning;   a  is  a  one -dimensional  array 
containing  10  elements  which  are   referred  to  by  a[l]^  a[2],  ...,  a[lO],  and  b  is  a. 
two-dimensional  array  containing  4x3  elements  which  are  referred  to  by  b[l, 3]> 
b[l,4],  b[l,  5];i  b[2, 3]>  "oo^  h[4, 5]»  In  an  ALGOL  program,  the  subscripts  may  be 
arbitrary  arithmetic  expressions  (See  below),  in  which  case  the  value  of  the 
expression  (if  necessary  properly  rounded  to  the  nearest  Integer)  determines  the 
particular  element  of  the  array. 

Nothing  is  said  in  ALGOL  about  the  fo:nn  in  which  Information  is  stored 
in  the  storage  locations.  Different  types  of  information  may  be  specified, 
however,  by  declarations  for  the  storage  locations.  Three  types  of  information 
are  possible;  real,  integer,  boolean.  A  declaration  real  a,,b  means  that  storage 
locations  a  and  b  contain  a  "real  number"  and  a  declaration  boolean  c,d  means 
that  storage  locations  c  and  d  may  contain  only  boolean  values  true  and  false. 
In  an  ALGOL  program,  a,  declaration  must  be  given  for  each  variable  and  each 
array  (in  the  latter  case,  the  type  is  put  in  front  of  the  array  declaration; 
in  case  real  the  type  may  be  omitted) , 


7.  This  means  "a  number  in  floating  point  representation". 
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Nothing  can  be  said  in  ALGOL  about  the  numerical  addresses  which  are 

associated  with  a  symholic  address,  and  there  is  no  direct  way  of  using  numerical 

8,9,10 
addresses.  '  ' 

The  main  constituents  of  ALGOL  programs  are  statements.  A  statement  is 
a  self  contained  command  to  do  something.  There  are  several  types  of  statements, 
of  which  the  following  three  are  in  a  certain,  sense  basic: 

1)  assignment  statements;, 

2)  go  to  -  statements,  r.nd 

3)  conditional  statements. 

The  following  statements  do  not  really  add  to  the  scope  of  the  lauiguage,  but  add 
to  the  convenience  in  its  usage.  These  ares 

k)  compound  statements, 

5)  loops, 

6)  blocks,  and 

7)  procedure  statements. 

An  ALGOL  program  is  a  sequence  of  statements,  separated  by^  and  preceded  by 
declarations  concerning  the  type  of  variables  and  the  type  and  size  of  the 
arrays  occurring  in  the  program. 

A  few  words  about  the  different  statements: 

1.  Assignment  Statements 

i-rithmetic  and  booleaji  calculations  can  be  performed  with  the 
numibers  in  storage  locations  and  with  constants,  and  this  is  indicated 
in  the  usual  way.  The  result  of  these  operatioiis  may  be  stored  and  this 
is  indicated  by  the  symbol  :  :=  meaning:  The  -^ralue  to  the  right  of  :  = 
is  stored  In  the  storage  location  to  the  left  of  :  =  .  Possible  assign- 
ment statements  8.re: 


8,  This  may,  however,  be  done  by  "procedures  written  in  machine  code". 
See  below. 

9o  It  is  not  assumed  that  the  same  numerical,  address  fee  associated  with  a 
symbolic  address  throughout  the  e.xec'ation  of  the  program.  On  the  con- 
trary, e.g.  in  blocks,  dynataic  storage  allot^ation  is  a"! most  necessary 
to  make  effective  use  of  the  storage  of  the  machine. 

10.  It  is  also  not  possible  to  state  whether  a  symbolic  address  refers  to 

storage  location  in  high  speed  memory  or  on  drums  or  tapes.  See,  however. 
Section  3ol. 
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=  X  +  1 

=  (a  +  b)  X  (c[l]  +  d[i,j]) 

=  a  <  0 


In  the  last  case,  b  must  be  declared  boolean  in  the  program  in  which  the 
statement  occurs,  and  the  assignment  statement  assigns  b  the  value  true 
or  false  depending  on  the  comparison  on  the  right  side  of  :  =  . 

X:  =  a  V  b  /s  c  +  d  <  e 
In  this  case,  a,  b  and  X  are  boolean,  and  c,  d,  e  may  be  integer  or  real. 
The  following  statements  are  not  assignment  statements? 

a+b;  =c+d(no  storage  location  on  left  side) 

Xt   =  X  <  y  (if  X  is  boolean,  it  cannot  occur 
in  a  comparison,  and  if  X  is 
arithmetic,  it  must  not  be  assigned 
a  boolean  value). 

2o  Go  to  Statements 

The  statements  of  an  ALGOL  program  are  executed  in  the  sequence 
in  which  they  are  written  down.  Some  of  the  statements  may  have  labels 
(which  have  the  form  of  identifiers  or  which  are  integers),  and  statements 
with  labels  again  may  have  labels.  The  label  is  separated  from  the  state- 
ment by  a  colon.  A  statement  -  go  to  L  -  where  L  is  some  label  means  that 
the  statement  with  label  L  is  the  next  one  to  be  executed. 

There  is  also  the  possibility  of  going  to  one  of  several  labels, 
depending  on  the  value  of  an  arithmetic  expression,  where  these  labels 
have  been  listed  and  given  a  name  in  a  switch  declaration. 

Example"   switch  alpha:  =  LI,  L2,  k,   Mj 
go  to  alpha  [k]^ 

The  go  to  statement  goes  to  statement  with  label  LI,  L2,  4,  M  if 
the  value  of  k  is  1,  2,  3?   ^  respectively. 


T  means  exponentiation.  The  symbol  I  i 

LS  e. 
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11.  The  arrow  I  means  exponentiation.  The  symbol  i  is  stronger  than 
a  I  b  +  c  means  a  +  c. 


+, ^ -,  X,  /  in  the  same  sense  as  e.g.  x  is  stronger  than  +.  Thus 


3.  Conditional  Statements 

These  make  the  execution  of  a  statement  (often  a  go  to  statement) 
dependent  on  the  validity  of  a  condition.  The  condition  is  a  boolean 
expression,  (See  assignment  statement). 

Condition  statements  have  one  of  two  forms s 

1)  if  condition  then  statement; 

2)  if  condition  then  statement  else  statement; 

Examples; 

if  a  <  b  then  go  to  L; 

ifa<bA  c  Ad  then  a;   =  b  else  Xs   =  Y; 

The  statements  depending  on  the  condition  may  not  themselves  be 
conditional  statements.  The  statement  after  the  else^in  forml)  may  be 
conditional  itself  in  either  of  the  two  forms. 

Example  s 

if  a  <  b  theft  &2-^  ^  else  if  a  =  b  then 
go  to  N  else  go  to  M; 

An  ALGOL  Programs 

integer  i;  array  a[lslO] 

is  =  1; 

L2  =  a[i];  =  i  |  2; 

is  =  i  +  Ij  if  i  >  10  then  go  to  L; 

k.      Compound  Statements 

Several  statements  may  be  combined  to  form  a  single  statement  by 
enclosing  them  with  the  statement  parenthesis  begin  and  end .  This  is  e.g. 
useful  when  a  sequence  of  statements  is  made  dependent  on  a  condition. 

Example i 

if  X  <  Y  then  begin  h:  =  X;  X;  =  Y;  Ys  =  X  end; 

This  little  program  assures  that  after  its  execution  the  number  in  storage 
location  X  is  greater  that  or  equal  to  the  number  in  storage  location  Y. 
Leaving  the  statement  parenthesis  begin  and  end  oat  would  make  only  the 
statement  h:  =  X  dependent  upon  the  condition.  The  little  program  could 
be  written  without  the  statement  parenthesis  in  the  following  ways 
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if  X  >  Y  then  go  to  L| 
h;   =  X|  X:   =  Yj   Ys   =  X} 

Lt      (Next  statement) 

There  are  other  uses  for  the  compo-und  statement,  e,g".  loops. 


Loopi 


There  are  three  types  of  loops,  which  we  call  list-type,  progress ion - 
type,  and  condition -type .  Generally  speaking,  a  loop  causes  a  statement 
(single  or  compound)  to  be  executed  several  times,  each  time  with  a  different 
value  of  the  loop  -variable .  In  the  list-type  loop,  the  values  of  the  loop 
variable  are  determined  by  a  list  of  arithmetic  expressions.  In  the  progres- 
sion-type loop,  the  values  of  the  loop  variable  are  determined  by  initial 
value,  step  width  and  final  value  of  an  arithmetic  progression.  In  the 
condition- type  loop,  the  values  of  the  loop  'mriable  are  determined  by 
evaluating  an  arithmetic  expression  on  each  pass  through  the  loop,  where 
this  expression  contains  the  loop  variable.  The  exit  from  the  loop  is 
made  when  a  certain  condition  is  fulfilled. 


Examples ; 


List -type  loops 

for  ei.1   =   loYj)  ^-S;.  9»6,  3=2  do  print  '"  a  T  2; 

Progression-type  loops 

for  is  =  1  step  1  until  10  do  a[i];  =  i  t  2; 

Condition-type  loop; 

X;  =  aj  d;  ==  0.5  x  (i-a); 
for  X;  ^  X  +  d  while  abs(d  t  2/x)  >  10   d£ 

ds  =  0.5  X  (a/x  +  X)^ 


The  la,st  statement  means  in  fulls 

Ls  Xs  -  X  +  d| 


if  abs(d  t  2/X)  >  10   then  begin 
ds  =  0.5  X  (a/x  +•  X)j  go  to  L  end; 

So  the  program  computes  -v  a  with  a  relative  accuracy  of  2  x  10 


12.  The  word  print  is  an  address  location  invention  (it  is  not  in  ALGOL  6o) 
with  the  meanings  Print  the  value  of  the  following  expression  or 
expressions  on  the  typewriter. 
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The  general  type  of  the  for  statements  can  be  easily  extracted 
from  these  examples.  It  might  be  worthwhile  to  add  that  there  are  no 
restrictions  on  the  kind  of  arithmetic  expressions  ,.^)oolea.n  expressions, 
or  statements  at  the  proper  pll.aces  In   loops.       '  « 

Also  mixed  loops  can  be  formed. 
Example  s 


This  is  equivalent  to  the  following  statement? 

k 

for  1;  =1  step  1  imtil  k-1^  k-1  step  2  imtil  10   do  S; 

This  in  turn  is  eqiii-valent  to  the  following; 


for  is  =  1  step  1  until  k-1  do  Sj 
for  i I   =  k-1  step  2  until  10  ^  S| 


Blocks 


A  block  is  a  compound  statement  that  contains  at  its  beginning 
declarations  about  some  of  the  quantities  used,  in  it.  Those  quantities  which 
are  declared  in  the  block  are  locg.l  to  that  block.  This  means; 

1)  They  are  not  accessible  in  statements  outside  the  block. 

2)  They  are  different  from  any  quantities  outside  the  block, 
even  if  the  same   identifier  is  used  for  quantities  outside 
the  block. 

Those  quantities  which  are  not  deelpir'^d  inside  the  block  are  called  global  to 
the  block. 

7.  Procedure  Statem.enbs 

These  are  best  discafesed  in  connection  with  procedure  declarations. 

Procedures  are  the  subroutines  of  ALGrC'L.  Those  ALGOL  programs  which 
are  believed  to  be  asefal  in  several  connections  with  maybe  different  "para- 
meters" each  time  may  be  declared  proced:;ireg  by  preceding  them  with  a  procedure 
heading.  In  this  heading,  a  name  i.s  given  to  that  program  and  the  parameters 
of  the  program  are  listed.  Farthemnsore,  certain  assertions  as  to  the  kind  and 
type  of  parameters  may  be  mad.e  in.  the  heading.  No  rigid  rules  for  the  pro- 
cedure heading  will  be  given  here. 
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In  a,  procedure  statement  (or  procedure  call)  the  name  of  the  procedure 
is  given  together  with  a  list  of  "actual  parameters".  This  call  is  equivalent 
to  execution  of  the  program  corresponding  to  the  name  mentioned  in  the  procedure 
statement,  after  suitably  replacing  the  "formal  parameters"  of  the  procedure 
declaration  by  the  actual  parameters  of  the  procedure  call. 

Further  discussion  should  be  given  to  parameters .  The  actual  para- 
meters may  be  substituted  for  the  formal  parameters  of  a  procedure  are  of  the 
following  type: 

1)  Values  which  are  computed  in  the  calling  program 

2)  Names  of  procedures  of  the  calling  program 

3)  Labels  of  the  calling  program 

k)     Storage  locations  of  the  calling  program 

5)  Rules  for  computing  a  value 

6)  Rules  for  computing  a  storage  location  (restricted  to 

subscripted  variables) 

7)  Names  of  arrays 

Furthermore,  if  a  procedure  has  a  single  result  like  a  procedure 
for  square  root  or  determinant,  this  result  can  be  omitted  from  the  parameters. 
Such  a  procedure  must  then  be  called  by  using  it  inside  an  arithmetic  expression, 
and  the  "result"  of  the  procedure  is  used  in  evaluating  that  expression.   In 
this  case,  the  type  of  the  result  must  be  given  in  the  declaration  of  procedure. 

¥e  give  some  examples  for  procedures  using  different  kinds  of  para- 
meters . 

Example  1  (parameters  of  type  l); 

real  procedure  sqrt  (a, n)j  value  a,nj  real  a;  integer  n; 
begin  sqrt  Is  =  aj  for  is  =  1  step  1  until  n  do  sqrt  1;  = 
0.5  X  (a/sqrt  1  +  sqrt  l)  end; 

The  parameters  a  and  n  are,  because  of  the  value  part  in  the  pro- 
cedure heading,  replaced  by  values  of  the  calling  program. 

Example  2  (parameters  of  type  2)1 

— —  ^b 

Integration  of  a  function  _yg   f(x)  dX  with  n  steps  using  the 
trapezoidal  rule* 
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real  procedigre  Int  (ajb^n.f)^  value  a,b^n^  real  a^b;  Integer  n; 
begin  integer  ij  ints  =  0; 

for  i;  =  1  step  1  until  n-1  do  ints  -   int  +  f(a  +  i  x  (b  -  a)/n); 
ints  =  (2  X  int  +  f(a)  +  f(b))  x  (b  -  a)/n  end; 

The  parameters  a^b  and  n  are  of  the  kind  already  discussed.  The 

parameter  f  is  a  procedure  with  one  V8.riable  as  input  a,nd  one  value  as  output. 

The  procedure  int  itself  ::a,s  one  single  result  as  output,  and  thus  is  called 
by  using  it  in  arithmetic  expressions. 

Example  3  (parameters  of  type  3^  '+)• 

procedure  sqrt  2   (a, n,X, L);  value  a,n|  real  a;  integer  n; 
if  a  <  0  then  go  to  L  else  Ss  s.  sqrt  (a,,n)j 

Parameters  a  and  n  are  of  the  typ^h;  already  discussed.  Parameter  X 
is  a  storage  location  of  the  calling  program  and  the  square  root  will  be  placed 
there  if  the  value  of  a  is  greater  than  zero.   If  a  <  0,  the  procedure  sqrt  2 
will  return  to  the  statement  with  label  L  (or  what  is  substituted  for  L  in  the 
call)  in  the  calling  program. 

Example  k   (parameters  of  type  5)s 

procedure  sum  (a,b,E,i)|  value  a,"bj  integer  a,bj 

begin  sum;  =  0;  for  is  -  a  step  1  uritil  b  do  sum:  =  sum  +  E  end; 

Parameters  a,b  (type  l)  and  1  (type  it)  are  already  discussed.  Para- 
meter E  is  of  type  5o  It  may  be  replaced  by  aa  arithmetic  expression.  While 
this  is  also  true  for  the  value  type  parameters^  in  the  latter  case  the  value 
of  the  expression  is  evaluated  before  entering  the  procedure.  While  in  the 
non-value  case,  the  rule  for  evaluating  the  expression  as  a  whole  is  entered 
as  actual  parameter  into  the  procedare. 

Examples s 

sijm  (l,10,a[j],j) 
sum  (l,n,  i  '[   2,  i) 
sum  (^,100.,  1  f  k  t  2,k) 
Another  example  could  be: 

sum  (l,  10,  sum(l,10,a[i,k]  j  2,i),k) 
which  would  evaluate  the  double  sum:         ^„ 

E   Z   a[i^k]  '  . 


In  this  case,  procedure  sum  calls  Itself  while  still  active.  It  must  he 
coded  accordingly.  These  "recursive  procedures"  (recursive  either  by- 
parameter  setting  or  by  fixed  procedure  calls)  are  as  convenient  for 
stating  computation  processes  as  they  are  bad  for  efficient  machine 
translation. 

Example  ^  (parameters  of  type  6}t 

Normalization  of  vectors  or  rows  or  columns  of  matrices, 

procedure  Norm  (a,n, l);  value  n;  Integer  n; 

begin  real  s; 

St   =  0|  for  is  =  1  step  1  -jntll  ndos;  -s  +  a|  2; 

for  is  =  1  step  1  until  n  do  as  =  a/sqrt  (s)   end; 

Examples  for  procedure  calls; 

Norm  (a[i],n,i) 

Norm  (a[j,k],n,j) 

Norm  (a[k,  j],n;,  j) 

Norm  (a[k,k],n, o) 

Norm  (a[k, j, j],n;, j)  (if  the  columns  of  the  matrix  a  are 

stored  in  some  strange  way.) 

Example  6  (parameters  of  kind  7)° 

Sum  of  elements  of  a  vector. 

procedure  sum  1  (a,n);  value  n;  integer  n| 

begin  integer  i  sxan  Is  =0; 

for  is  =  1  step  1  until  n  do  s-jw.  l.t  ■-■  s'jm  1  +  a[i]; 

This  procedure  is  more  restricted  than  procedure  sum,  but 
probably  faster. 


13.   If  the  translator  is  not  general  eaough  to  find  out  that  s  does  not 

change  in  tha^  loop,  and  codes  accordingly,  1~   would  be  wise  to  rewrite 
it  as;   s;  =  sqrt  (s);  for  1:  =  1  step  1  until  n  do  as  =  a/s. 
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,2«2  Sorting  of  the  Rows  of  a  Jfetrix  According  to  the  Size  of  the  First  Element 
of  a  Row 

METHOD; 

Given  a  matrix  A  of  order  m,  n  with  elements  a(i,  j),  this  program 
jlnterchanges  the  row  so  the  a(i,l)  >  a(i+l,l)  for  i  =  1,2, ,.,m-l. 

Rather  than  work  with  the  complete  rows  in  finding  the  proper  order, 
we  first  determine  the  desired  order  of  the  rows  by  generating  two  arrays, 
i[j]  =  f  and  k[p]  =  q,  in  three  steps. 

i[j]  =  f  means  row  member  f  goes  into  row  number  j. 

^[p]  =  ^  means  row  nxanber  p  goes  into  row  number  q. 

STEP  1:  Array  i[j]  starts  out  with;   i[l]  =  1,  i[2]  =  2, ...i[m]  =  m. 

STEP  2;  Consider  j  =  1,2,  ...,m-l  until:  a(i[  j+l],l)>a(i[  j],l) . 

Interchange  the  value  of  i[j]  in  zhe   array  with  the  value  of 
i[o+l]«  Subtract  one  from  last  j  unless  j  =  1  and  test  the 
above  inequality.   If  it  holds, make  the  interchange  of  values 
and  continue  with  the  subtracting  and  testing.   If  it  doesn't 
hold, or  if  j  -  1,  repeat  the  process  of  considering  j  =  1,2, ,..,m-l 
until  the  inequality  holds  again  or  until  j  =  m-1.  Now  the 
i[j]  =  f  array  is  completed,  , 

STEP  3 2  The  k[p]  =  q  array  is  formed  by  setting  k[i[j]]  =  j  for 
J  =  1,  (d,,«,,m, 
e  ,g. 


STEP  3: 


3 

i[l] 
i[2] 
i[3] 
±[h] 
i[5] 


1  2 


=  d.       d 

—  X 

=  3 

=  k 

=  5 


2 

3 
If 
1 
5 
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Once  the  order  is  determined  we  rearrange  the  rows  in  this  manner^ 

a)  The  first  row  that  needs  interchanging  is  put  in  temporary  storage. 

b)  From  temporary  storage  this  row  goes  to  its  proper  row  number  and  the 
row  that  is  replaced  is  put  in  the  aforementioned  temporary  storage. 
(Change  the  k[  ]  array  for  the  row  number  which  is  in  its  proper  position 
to  be  equal  to  -1.) 

c)  Repeat  process  b)  for  the  row  now  in  temporary  storage  until  the  end  of 
a  cycle  is  reached,   (e.g.:  row  #1—^  row  #9  — ^  row  #5  -->row  #l) 

d)  At  the  end  of  a  cycle  the  row  in  temporary  storage  is  put  in  its  proper 
position  without  the  row  being  replaced  put  in  temporary  storage  since 
it  has  already  been  repositioned. 

e)  The  "new"  matrix  is  now  checked  for  the  first  row  out  of  position  and 
the  process  a)  through  d)  is  repeated.  If  none  of  the  rows  are  out  of 
position  {e.g. I     k[l]  =  k[2]  =  ...  =  k[m]  =  -l)  stop. 

RI  is  the  number  of  the  row  being  moved  to  its  proper  position. 

RO  is  the  number  of  the  row  replaced  by  row  number  RI. 

Q[l]^  Q[2], . . .,Q[n]  is  used  as  temporary  storage  locations  for  a  row  of  the 

matrix  A. 
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FLOW  CHAET 


START 


->- 


Form  i[  ]  array; 


Order  the  i[  ]  array; 


V 


Form  the  k[  ]  array; 
k[p]  =  q. 


Start  with  the  row 
number  equal  to  Ic 
(Rl  =  1) 


Calculate  which  row 
position  row  number 
RI  goes  into. 


i  Preset  test  for  the 
I  completion  of  row 
.    number  RI. 


Is  row  number  RI 
already  in  the  right 
position? 


No 


^/   Yes 


Take  the  first  row  in 
the  new  matrix  which 
is  out  of  position 
and  assign  its  row 
number  to  RI. 

I\ 


STOP 


< 1 

/,   Yes 

Are  all  of  the  rows  in 
the  correct  position? 

No 

^ 



Put  the  row  in  temporary 
storage  in  its  correct 
row  number.   (RI-^RO) 


/\   Yes 


Does  row  number  RI  in 
temporary  storage  go  into 
a  row  location  that  has 
already  been  completed? 


No 


^- 


/\ 


Preset  test  for  the 
completion  of  row 
number  RI. 


Renumber  rows  making  the 
new  RI  equal  the  old  RO. 


Put  the  row  in  temporary 
storage  in  row  position 
nijmber  RO  and  put  row 
number  RO  in  temporary 
storage . 


/s 


Put  row  number  RI  in 
temporary  storage 
Q[l]  =  a[RI,l], 
Q[2]  =  a[RI,2],..., 
Q[n]  =  a[RI,n]. 


-« 
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ANNOTATION 


for  j ;  =  1  step  1  until  n  do  i[j]:  =  j; 


STEP  1,  See  METHOD. 


for  j;  =  1  step  1  imtll  n-1  do  if  a,[i[j+l],l]  >  a[i[j],l] 

then  for  kt   =  j+l,  k-1  while  a[i[k],l]  >  a[i[k-l],l] 
do  begin  h;  =  i[k3i  i[k];  =  i[k-l];  i[k-l]:  =  h  end; 

for  j2  =  1  step  1  until  mdok[i[j]]o  =j; 


STEP  2,  See  METHOD. 


STEP  3>  See  METHOD. 


Ell   =  1; 


St&rt  by  putting  row  #1 
in  the  right  order. 


LI:  RO;  =  k[Rl]| 


Row  #RI  goes  into  row  #R0. 


k[Rl];  =  -Ij 


Preset  test  for  completion 
of  row  #RI. 


if  RO  =:  RI  then  go  to  I2j 

for  js  =  1  step  1  until  n  do  Q[j]:  =  a[RI,j]j 


Row  #RI  need  not  he  movedo 

Row  #RI  needs  to  he  moved 
so  put  row  #RI  in  temporary 
storage . 


L°  fo^  js  =  1  step  1  until  n  do 

begin  hi   =  a[RO,j]j  a[RO,j];  =  Q[j]i  Q[j];  =  h  er^ 

Rli   =  RO;  ROs  =  k[Rl]; 
k[Rl]:  =  Ij 


Interchange  the  corresponding 
elements  of  temporary  storage 
and  row  #R0. 


Consider  a  new  RI  and  RO. 


Preset  test  for  completion 
of  row  #RIo 


i 


if  k[RO]  =  -1  then  for  j  s  =  1  step  1  until  n  do 
begin  a[RO,o]:  =  Q[j3i  go  to  L2  end 


If  at  the  end  of  a  cycle 
put  the  temporary  storage 
row  in  the  right  row 
location. 


else  go  to  Lj 


If  cycle  isn't  completed, 
repeat . 


12:  for  j:  =  1  step  1  until  n  do  if  k[j]  ^  -1  then 
begin  RI:  =  jj  go  to  LI  end 


else  stop. 
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Set.  RI  equal  to  the  first 

row  of  the  matrix  which  is     J 

now  out  of  order  and  repeat.   S 


None  of  the  rows  are  out 
of  order  so  stop. 


i 


c.tj      a^jcii  DC  1'. 


la,  oj.  XI- CO 


STORING  A  SPARSE  MATRIX 


METHOD; 

A  matrix,  A,  of  order  m,n  is  read  ir_  by  a  seq\;.ence  of  triplets  - 
i,j,a(l,j)  -  where  a(i,j)  ^  0,  The  last  triplet  Is  followed  by  0,0,0  to 
indicate  the  end  of  the  triplets,  k  ccxits  the  number  of  non-zero  elements 
of  A  read  in.  Three  arrays  are  generated  -  XI [k],  X2[k],  X3[i].  Where 
a(i;j)  is  the  k"  element  read  in,  Xl[k]  -  a(i,  j),  X2[k]  =  j;  and,  if  a(i,j) 
is  the  first  non-zero  element  of  a  row,  X3[i]  =  k.  If  a  row  contains  no 
non-zero  element,  X3[i]  =  0. 


Thus,  the  element  in  XI [k]  lies  in  col^jmn  number  X2[k]  and  in  row 


i  where: 


X3[i]  <  k  <  X3[i/l/q],  q  -  0,1, ...,m-l  until  X3[i/l/q]  /  0. 

For  this  reason  we  set  X3im/l]  ~  total  number  of  non-zero  elements  /I. 
I  is  the  number  of  the  row  being  worked  upon. 


FLOW  CHART 


START- 


Read  in 
first 

triplet  < 


Set: 
X3[I] 


Increase 
row 

n.-'iniber; 
(l-~l/i). 


1^   Yes 


Test  fo2 
new  row. 


i 


No 


Start  count 
I  of  non-zero 
I  el.ements;  k=l 

^ 


No 


> 


Test  if 

a(i,j)  is  in 
row  I,  (  L~l) , 


\^ 


Set:  XIFIJ  -  kl 


y 


Set : 

Xl[k]  -  a(i,j); 

Xt[k]  -  j. 


1 


No  : 


Increase  numbe. 
of  non-zero 
elements  read 
in,  (k  =  k/l). 


N^ 


Test  foic 
end  of 

triplets . 


Read  in  new 

triplet „ 


t 


Yes 


Start  count 

of  row 
n  amt  er;  I--1 . 


Increase 
row  number, 
j  (I  =  l/l). 


->- 


No 


No 
> 


Set:  X3[I]=C 
for  new  row 
number . 


> 


STOP' 


^ 


Test  for  zero 
matrix. 


^'  Yes 


Set:  X3[I]  =  0, 
where  1=1. 


\' 


Test  for  m 
row,  (l  =  m) . 


sf    Yes 


< 


Set  dummy  row  array 
element:   X3[m/l]=k, 
where  k  at  this  point 
equal  the  total  number 
of  non-zero  elements 
plus  one. 
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PROGRAM' 


Ik 


ANNOTATION 


Read   i,  j,a; 
ks  =  Ij 
I:  =  1; 

zero:  if  i  =  0  then  "begin 

X3[l]s  -  Oi 
go  to  fin  end; 

test;  if  i  j^  I  then  "begin 
X3[l]o  -  0; 
go  to  adjust  endj 

X3[l]s  -  kj 

go  tc  entry; 


adjust; 


I:  =  I  /  Ij 
go  to  test; 


entry?  Xl[k];  =  a; 
X2[k]s  =  j; 
k?  =  k  /  1; 

if  i  =  0  then 
go  to  fin; 

if  i  ^  I  then 
go  to  adjust; 

go  to  entry I 


fin;  if  I  ==  m  then  begin 
X3[l/l]s  -  k; 
16 


stop" 


end.; 


Is  =  I  /  1; 
X3[l]s  =  0; 
go  to  fin; 


First  triplet. 
Count  first  element. 
Start  with  first  row. 

This  handles  the  case  of 
the  zero  matrix. 


This  takes  care  of  the  rows 
which  contain  no  non-zero 
elements . 

a( i, j )  is  the  first  non-zero 
element  of  the  row. 

The  number  of  the  row  "being 
worked  upon  is  increased  "by  1. 

Entry  made  to  element  array. 
Entry  made  to  column  array. 
Increase  coiait  on  non-zero  elements, 

Read  in  new  triplet. 

Test  for  the  end  of  the  triplets. 

Test  for  the  start  of  a  new  row. 


a(l^, j)  isn't  the  first  element 
of  a  row. 

Dummy  row  array  element  is  set 
up  so  that  t'he  end  of  row  m  can 
be  determined. 

This  completes  the  row  array 
if  the  last  rows  contain  no 
non-zero  elements. 


1^.  No  declarastiions  are  given. 

15.  Read  is  a  symbol  added  tc  ALGOL  for  the  purposes  of  this  report.  It  has 
the  meaning:  Read  the  numbers  (ALGOL-Form)  on  standard  input  mediTan  and 
store  them  in  the  storage  locations  given  in  the  read  statement, 

16.  Stop  is  a  symbol  added  to  ALGOL  for  the  purpose  of  this  report;  meaning  is 
obvious . 
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MULTIPLICATION  OF  TWO  SPARSE  MATRICES 

METHOD; 

Two  spaxse  matrices,  X  and  Y^,  are  stored  aacording  to  the  program, 
"Storing  a  Sparse  Jfetrix".  X  is  of  order  m,n  and  Y  is  oi"  order  n,p.  The 
three  arrays  for  matrix  X  are  XI  [  ],  X2[  ],  and  X3[  ]»  The  three  arrays  for 
matrix  Y  are  Yl[  ],  Y2[  3,  and  Y3[  ]. 

The  first  column  of  Y  is  obtained  usirig  the  three  triplets  of  Y  and 

is  stored  in  temporary  storage;  Y(  1  ),■.<,  o,Y(n)o 

The  first  row  of  X  Is  multiplied  Iby  the  col-umn  of  Y  in  temporary 
storage.  If  the  element  of  Z,  Z(i,j),  thus  formed,  is  the  f   non-zero  element 
of  Z,  Z(i, j)  is  stored  hy  a  triplet  Zl[fl,  Z2.U,},   Z3[f],  where  Zl.[f]  =  Z(i,j), 

Z3[fl  =  i|  aad  Z2[j]  ^   f  if  Z(i, j)  is  the  first  aon-zero  element  of  a  columno 

After  all  of  the  rows  of  X  haye  been  m-uJ.tipl:led  by  the  column  of  Y 
in  temporally  storage,  the  next  c^lijann  of  Y  is  stored  in  the  temporary  storage 

and  the  mialtiplication  starts  overo 

CY  is  the  niMber  of  the  column  being  worked  upon  of  Y. 

RY  is  the  number  of  the  row  being  worked  upon  of  Y. 

RX  is  the  number  of  the  row  being  worked  upon  of  X. 

f  is  used  to  count  the  number  of  .uoa-zero  elements  of  Z 
as  they  are  formed » 

E  is  used  as  temporary  sto;rage  .in  fo;ncliig  an  elonent  of  Z. 
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FLOW  CHART 


Set.  number  of 
non-zero  ele- 
ments of  Z  to  If 


{■■ 


1 


\/ 


START 


STOP 


Set  the  dummy 
element  of  the 
column  array  of 
Zi      Z2(p+l)=f, 
where  f  is  equal 
to  the  total 
nomoer  of  non- 
zero elements 
cf  Z  plus  oneo 


Yes 


Is  this  not  a 
last  column  of 
Y,  (CY  =  p+l)? 


Ko 


Set  the  column 
nimiber  of  Y 
equal  to  1, 
(CY  =  1). 


<■ 


Yes 


Set  h  =  Ij 
(after  first 
element  of 
column  is 
stored,  h  is 
set  equal  to 
zero ) . 


Set  next  ele- 
m.ent  of  column 
array  of  Z 
equal  to  zero. 


Set  row  number 
of  X  equal  to 
1,  (RX  =  l). 


N/ 


Clear  sum  lo- 
cation for 
element  of  Z, 


No 


Is  this  not  a 
last  row  of 
X,  (RX  =  m+1)? 


v^ 


Yes 


Increase 
column  number 
of  Y  (CY  =  CY+l). 


Set  the  row 

number  of  Y 

equal  to  1, 

(RY  =  1). 


/\ 


/\ 


No 


Is  this  not  a 

row  of  Y, 
(RY  -  n+1)? 


Nc 


Increase  row 
nimiber  of  Y, 
(RY  =  RY+l). 


<■ 


-<- 


Is  Z(l,^)  -  0? 


No 


Multiply  the 
column  of  Y  in 
temporary  stor- 
age by  a  row  of 
X^  keeping  the 
simi  in  Ea 


A  No 


Is  present  row 
of  X  all  zero 
elements? 


>/  Yes 


Increase  row 
niaflber  of  X, 
(RX  =  RX+l). 


-^ 


Yes 


No 


Set  next  temporary 
storage  location 
equal  to  zero, 
(Y[RY]  =0). 


T 


Do  arxy  of  the 

elements  in  row 

RY  lie  in  column  CY? 


Yes 


Store  that  element 
in  Y[RY]. 


Set: 

Zl[f]  =  E, 
Z3[f]  =  RX. 


Is  Z(i,j)  the  first 
non-zero  element  in 
the  column  of  Z, 
(h  -  1)? 


^'Yei 


Set: 


Z2[CY]  =  f, 

h  =  Co 


Increase  the  covint 
of  non-zero  elements 
of  Z,  (f  =  f+l). 
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ANNOTATION 


for  CYi   =   1  step  1  i;uitll  p  do 
begin  for  RYs  =  1  step  1  imtll  n  do 

begin  Y[RY];  =  0;  for  b;  =  Y3[RY3 

step  1  imtll  Y3[RY+1]  -1  ^ 
begin  if  Y2[b]  =  GY  then 

begin  Y[RY]s  =  Yl[b]j  go  to  T  end 
else  if  Y2[b]  >  CY  tjien  go  to  T  eadj 


Set  count  of  non-zero 
elements , 


Stores  one  colxann  of  the 
matrix  Y  in  temporary 
storage , 


Ts  endj 
h:  =  Ij 

ze[cY]s  =  0^ 

for  RX;  =  1  step  1  mi.tll  m 
begin  Es  ~  Oj 

if  X3[RX]  ^   0  then 

begin  for  k^  ^^   X3[RX]  step  1  ■until  X3[RX  +l]  - 
Es  r^   E  +  Xl[k]  X  Y[X2[k]; 
if  E  ~  0  then  go  to  D  else  go  to  C  end 
else  go  to  D| 


Preset  test  for  entry  to 
the  coliann  array,  and  preset 
the  column  array  element. 


Form  one  element  of  Z. 
do 


Co  Zl[f]s  -:  Ej  Z3[f]s  -  RX| 


if  h  =  1  then  begin  Z2[CY]:-fj  h^=  0  ends 


f  s  =  f  +  1; 


Entry  is  made  to  the  element 

array  and  the  row  array. 

Check  to  see  if  an  entry 

is  made  to  the  colimm  array. 

Count  non-zero  element. 


Dr  endj 
go  to  W| 

W;  endj 
Z2Ep+l];  . 
stop. 


Set  dummy  element  of 
column  array  so  that  the 
end  of  column  p  can  be 
determined  and  stopo 
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There  is  another  way  of  storing  the  information  contained  in  the  array  i. 
Essentially,  array  i  stores  sequences  of  increasing  integers,  A  good  way  to  do 
this  on  a  binary  ma,chine  would  "be  this;  a  two-dimensional  boolean  array  c  is 
formed,  an  c[k,  e]  is  true  if  the  j ""  element  of  row  k  is  not  zero.  Otherwise, 
c[k, e]  is  false .  The  boolean  array  is  stored  in  the  machine  by  rows,  with  one 
word  for  each  row. 

This  arrangement  can  be  described  in  ALGOL,  and  the  multiplication  of 
two  matrices  stored  in  this  form  can.  be  programmed  in  ALGOL.  No  attempt  is  made 
here  to  do  this.  The  question  of  whether  this  is  good  or  not  depends  on  the  trans- 
lator. If  it  is  clever  enough  to  use  individual  bits  of  a  word  to  store  boolean 
arrays,  it  might  not  be  too  bad.  Otherwise,  it  is  out  of  the  question. 

Also  in  the  first  case  it  is  doubtful  whether  a  program  resulting  from 
machine  translation  would  be  good.   It  would  in  many  places  refer  to  the  elements 
a[i, j].  Computation  of  addresses  in  this  case  consists  of  computing  the 
address  of  the  i   word  of  the  array  and  shifting  this  j  places.   "Recursive" 
address -computation  would  pick  out  the  elements  of  a[i, j]  one  by  one,  shifting 
one  place  at  a  time  -  and  this  again  is  out  of  range  for  reasonable  translators. 
Even  worse,  some  machines  might  have  "normalizing"  orders  which  allow  one  to  find  the 
next  non-zero  bit  of  a  word  thus  increasing  speed  considerably,  A  way  to  write 
an  efficient  AIGOL  program  using  these  features  is  the  incorporation  of  "procedures" 
written  in  machine  code.  The  following  procedures  would  serve  very  well  in  the 
previous  example; 

procedure  setbit  (a, i,E)j  value  E, 1^  Integer  i;  boolean  E; 
begin  This  procedure  sets  bit  number  i  in  storage  location  a 
equal  to  1  if  E  is  true  and  zero  if  E  is  false  end; 

boolean  procedure  bit  (a, i)|  value  i;  integer  i; 

begin  This  procedure  is  a  function,  the  value  of  which  is  true 

if  the  i   bit  in  storage  location  a  is  1  and  which  is 

zero  otherwise  end; 

In  the  example,  instead  of  giving  machine  code  for  the  procedure,  a 
verbal  description  is  used. 
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III.  LIMITATIONS  OF  ALGOL 


ALGOL  was  designed  for  the  description  of  numerical  processes,  and 
ail  of  its  features  and  constructs  were  created  to  make  description  of  numerical 
processes  convenient  and  machine  independent »  Some  price  had  to  be  paid  for 
"machine  independence"  (this  was  discussed  earlier)  and  some  price  had  to  be 

paid  for  the  conveniences,  the  latter  being  that  some  computing  processes 

17 
cannot  be  described  at  all  in  AKJOL.  The  basic  deficiencies   of  ALGOL  seem 

to  be  theses 

a)  Handling  of  information  which  is  non-nimerical 

b)  Reference  to  different  types  of  memory  arranged  in  hierarchy 

c)  Explicit  computation  of  numerical  addresses 

d)  No  possibility  to  state  possible  parallel  computations 

More  accidental  than  basic  is  the  lack  of  input-output  statements  in 
ALGOL o  Input -output  statements  could  be  easily  incorporated  in  the  same  "machine 
independent  level"  as  the  other  stsAements  of  ALGOL. 

In  the  following  serjtions,  items  a)  and  b)  are  examined  in  more  detail 
and  some  means  of  overcoming  the  difficulties  are  disc-ussed. 

3<.l  Handling  of  Non-Numerical  Infoinnation 

The  imits  of  information  which  ALGOL  recognizes  are  numbers  and 
boolesm  values.  There  are  ittany  processes  where  the  units  of  information  are 
quite  different  entitles.  They  may  be  mathematical  formulas,  or  positions  of 
a  chess  board,  or  properties  of  names,  or  names  or  arrangements  of  matrices  with 
many  zero  elements.  Where  this  originally  non-numerical  information  can  be 
described  conveniently  in  numerical  fo:i°m  (as  e.g.  in  the  case  of  sparse 
matrices.  See  2.3)  ALGOL  is  good  and  useftil.  There  are  cases  where  this 
(always  possible)  numerical  description  is  Inadequate,  and  for  two  reasons: 

1)  The  description  of  the  proceas  in  these  numerical 
te^rms  is  inconvenierit . 

2)  The  ALGOL  program  caanot  be  translated  into  an 
efficient  machine  program. 


17.  It  is  not  altogether  clear  whether  these  deficiencies  should  be  removed 
at  all,  that  is,  whether  "deficiency"  is  the  right  word.  Nobody  would 
say  that  an  automobile  that  cannot  fly  or  cannot  climb  to  the  top  of  the 
Nanga  Parbat  is  "deficient". 
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1  A 

If  the  information  is  e.g.  the  property  of  an  individual   with  values 
"single,  married,  divorced",  one  can  store  this  informa,tion  in  a  symbolic  address 
"family  status"  and  store  single  as  the  number  1,  married  as  nimiber  2  and  divorced 
as  number  3-  An  ALGOL  program  with  this  and  several  more  conventions  of  similar 
caliber,  containing  statements  of  the  fozm  "if  family  status  =  1  then  go  to  tax; " 
is  however  far  from  readable.  On  the  other  hand,  such  a  program  cannot  be  trans- 
lated into  an  efficient  machine  code  since  the  values  of  the  variable,  family 
status,  need  only  2  bits  for  encoding,  and  there  is  no  way  of  telling  this  to  the 
translator;  and  of  course,  no  "reasonable"  translator  will  be  able  to  find  this 
fact  by  itself. 

Similar  remarks  apply  to  other  kinds  of  essentially  non-numerical  and 
non -boolean  information. 

3.2  Reference  to  Different  Types  of  Stora,ge 

No  "reasonable"  translator  will  be  able  to  find  out  which  symbolic 
address  in  an  ALGOL  program  to  put  into  auxiliary  memory,  (A.M.),  and  which  to 
put  into  high-speed  memory,  (HSS),  in  an  efficient  way,  A  translator  may  very 
well  be  able  to  generate  the  sequence  of  instructions  necessary  to  transfer 
information  from  A.M.  to  HSS  once  it  knows  that  they  have  to  be  constructed. 
There  is  no  direct  way  in  ALGOL  to  tell  this  to  the  translator. 

A  simple  way  out  is  the  use  of  procedures  written  in  machine  language. 
We  will  illustrate  this  by  an  example  for  Illiac. 

We  assume  that  we  have  machine  code  procedures  readdrum  (a,b)  and 
writedrum  (a,b)  where  a  is  a  symbolic  address  and  b  is  an  arithmetic  expression. 
These  procedures  do  the  following s   For  readdrum  the  number  stored  at  drum 
address  b  (numerical  address  l)  is  transferred  into  stora,ge  position  a,  (symbolic 
address;  either  simple  address  or  element  of  an  array),  and  for  writedrum  the 
opposite  transfer  t.a,kes  place. 

We  want  to  store  a  symmetric  matrix  on  the  drum.   In  order  to  save 
space  we  do  not  store  the  lower  triangle  of  the  matrix.  We  want  to  be  able  to 

read  and  write  full  columns  and  rows  of  the  matrix  in  minimum  access  time.  This 

19 
can  be  accomplished  in  the  following  way;  '^ 


l8.  The  individual  itself  is  of  course,  represented  by  an  integer  number,  a 
common  practice  before  the  advent  of  computing  machinery. 

19 «  This  arrangement  is  used  in  several  codes  written  by  Dr.  Dickman  for  the 
Illiac . 
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A  sequence  of  locations  on  the  drum  can  be  called  in  minimum  access 
time  if  their  (numerical)  addresses  differ  by  multiples  of  6k.     The  arrow  in 
the  preceding  picture  shows  the  location  of  the  third  row  (or  column)  of  the 
matrix.  The  formula  for  computing  the  address  of  element  a[i,j]  is: 


loc  (a[i,j])  = 


loc  (a[l,l])  +  (j-1)  X  65  +  i-1,   if  i  <  j. 
loc  (a[l,l])  +  (i-1)  X  65  +  j-1,   if  i  >  j. 

If  the  number  of  rows  of  the  matrix  is  greater  than  Sk,   the  arrangement  has 
to  be  changed.  We  do  not  go  into  this  here. 

The  computing  speed  of  the  Illiac  is  not  high  enough  to  allow  computation 
of  the  next  address  daring  the  reading  or  writing  of  two  successive  elements  of 
a  row.  Therefore,  the  addresses  of  elements  of  a  row  or  coliomn  are  best  pre- 
computed  before  the  transfer  is  entered.  The  program  for  reading  the  1   row 
(or  column)  of  the  matrix  into  the  array  reads  as  follows: 

for  j:  =  1  step  1  until  i  do  a.d[j]:  =  A  +  (i-l)  x  65  +  j-1; 
for  jt   =   i+1  step  1  until  n  do  a.d[j]:  =  A  +  (j-l)  x  65  +  i-l; 
for  j:  =  1  step  1  \mtil  n  do  Readdrum  (d[j],  ad[j]); 

Here  A  is  the  storage  location  that  contains  loc  (a[l, l]),  and  ad[j]  is  an 

integer  array  for  intermediate  storage  of  the  addresses  of  the  elements  of 
,th 


the  i 


70W. 
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In  the  first  two  loops  the  computation  of  ad[j]  should  he  done  recursively. 
If  the  translator  does  not  do  that,  the  programmer  should  rewrite  the  first  loop  as° 

ad[l]:  =  A  +  (i-l)  x  65  -1; 

for  j:  =  2  step  1  -until  i  do  ad[j]?  =  dd[j-l]  +1; 

Using  symbolic  addresses  a,lso  on  the  dr-om  would  be  a  slight  "improvement" 
over  the  previous  scheme,  at  least  as  far  as  the  formulation  of  computation  processes 
is  concerned »  For  the  final  machine  program,  the  programmer  might  want  to  have 
control  over  the  actual  physical  arrayment  of  his  information  on  the  drum  because 
this  may  determine  the  computing  speed  of  the  process. 

The  best  wa,y  of  describing  usage  of  auxiliary  store  would  be  as  follows 
(but  this  is,  of  course,  beyond  the  scope  of  ALGOL).  Some  declarations  are  added 

to  the  language  concerning  the  way  in  which  the  information  is  to  be  stored  in  the 

20 
auxiliary  memories  of  the  computer.    At  the  beginning  of  each  block,  declarations 

are  given  concerning  the  particular  information  which  is  needed  in  HSS  in  that 

block.  The  translator  when  scanning  these  declarations,  generates  the  instructions 

for  transfer  of  information  to  HSS  at  the  beginning  of  each  block,  governed  by  the 

declarations  concerning  this  information.  At  the  end  of  a,  block,  information  is 

transferred  back  to  A.M.,  if  necessary.  This  scheme  would  have  the  following 

advantages; 

1)  It  is  mechanically  translatable. 

2)  It  gives  the  programmer  control  over  arrangement  of 
his  information,  but  relieves  him  from  the  detailed 
coding  of  the  transfer  of  Information. 

3)  A  program,  from  which  all  the  declarations  concerning 
auxiliary  storage  are  taken  away,  is  still  a  valid 
description  of  the  computing  process, 

k)     Programs  cannot,  of  course,  be  exchanged  between 
installations  with  different  types  of  auxiliary 
memory,  but,  by  changing  the  declarations  concerning 
aiixiliary  memory  one  gets  a  program  that  works  for 
the  machine  with  a  different  memory. 


20.  Such  a  declaration  could  e,g.  look  like  this:  real  drum  array  a[l; 100,1; 100], 
but  other  forms  are  also  possible.  No  specific  proposal  is  made  here. 


3.3  Numerical  Addresses 

In  some  cases  it  might  be  worthwhile  to  refer  to  numerical  addresses. 
No  possibilities  to  do  so  exist  in  ALGOL.  The  incorporation  of  symbols  like 
loc  (...)  meaning  "the  numerical  address  of  the  symbolic  address  in  the 
parenthesis"  and  cont  (...)  meaning  "the  number  in  the  memory  position  the 
numerical  address  of  which  is  given  by  the  arithmetic  expression  in  the 
parenthesis"  is  dangerous.  Some  translator  might  have  dynamic  storage  allo- 
cations and  it  is  difficult  to  see  how  programs  can  be  exchanged  which  use  loc 
and  cont  without  a  whole  set  of  rules  concerning  storage  allocation. 
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